Biomechanical Effects of Ti-Base Abutment Height on the Dental Implant System: A Finite Element Analysis.

This study aims to analyse, using a finite element analysis, the effects of Ti-base abutment height on the distribution and magnitude of transferred load and the resulting bone microstrain in the bone-implant system. A three-dimensional bone model of the mandibular premolar section was created with an implant placed in a juxta-osseous position. Three prosthetic models were designed: a 1 mm-high titanium-base (Ti-base) abutment with an 8 mm-high cemented monolithic zirconia crown was designed for model A, a 2 mm-high Ti-base abutment with a 7 mm-high crown for model B, and a 3 mm-high abutment with a 6 mm-high crown for model C. A static load of 150 N was applied to the central fossa at a six-degree angle with respect to the axial axis of the implant to evaluate the magnitude and distribution of load transfer and microstrain. The results showed a trend towards a direct linear association between the increase in the height of the Ti-base abutments and the increase in the transferred stress and the resulting microstrain to both the prosthetic elements and the bone/implant system. An increase in transferred stress and deformation of all elements of the system, within physiological ranges, was observed as the size of the Ti-base abutment increased.

Understanding the influence of surgical parameters on craniofacial surgery outcomes: a computational study.

Sagittal craniosynostosis (SC) is a congenital condition whereby the newborn skull develops abnormally owing to the premature ossification of the sagittal suture. Spring-assisted cranioplasty (SAC) is a minimally invasive surgical technique to treat SC, where metallic distractors are used to reshape the newborn's head. Although safe and effective, SAC outcomes remain uncertain owing to the limited understanding of skull-distractor interaction and the limited information provided by the analysis of single surgical cases. In this work, an SC population-averaged skull model was created and used to simulate spring insertion by means of the finite-element analysis using a previously developed modelling framework. Surgical parameters were varied to assess the effect of osteotomy and spring positioning, as well as distractor combinations, on the final skull dimensions. Simulation trends were compared with retrospective measurements from clinical imaging (X-ray and three-dimensional photogrammetry scans). It was found that the on-table post-implantation head shape change is more sensitive to spring stiffness than to the other surgical parameters. However, the overall end-of-treatment head shape is more sensitive to spring positioning and osteotomy size parameters. The results of this work suggest that SAC surgical planning should be performed in view of long-term results, rather than immediate on-table reshaping outcomes.

BaHf0.05Ti0.95O3 Ceramics from Sol-Gel and Solid-State Processes: Application to the Modelling of Piezoelectric Energy Harvesters.

A typical piezoelectric energy harvester is a bimorph cantilever with two layers of piezoelectric material on both sides of a flexible substrate. Piezoelectric layers of lead-based materials, typically lead zirconate titanate, have been mainly used due to their outstanding piezoelectric properties. However, due to lead toxicity and environmental problems, there is a need to replace them with environmentally benign materials. Here, our main efforts were focused on the preparation of hafnium-doped barium titanate (BaHfxTi1-xO3; BHT) sol-gel materials. The original process developed makes it possible to obtain a highly concentrated sol without strong organic complexing agents. Sol aging and concentration can be controlled to obtain a time-stable sol for a few months at room temperature, with desired viscosity and colloidal sizes. Densified bulk materials obtained from this optimized sol are compared with a solid-state synthesis, and both show good electromechanical properties: their thickness coupling factor kt values are around 53% and 47%, respectively, and their converse piezoelectric coefficient d33∗ values are around 420 and 330 pm/V, respectively. According to the electromechanical properties, the theoretical behavior in a bimorph configuration can be simulated to predict the resonance and anti-resonance frequencies and the corresponding output power values to help to design the final device. In the present case, the bimorph configuration based on BHT sol-gel material is designed to harvest ambient vibrations at low frequency (<200 Hz). It gives a maximum normalized volumetric power density of 0.03 µW/mm3/Hz/g2 at 154 Hz under an acceleration of 0.05 m/s2.

The Use of Virtual Tissue Constructs That Include Morphological Variability to Assess the Potential of Electrical Impedance Spectroscopy to Differentiate between Thyroid and Parathyroid Tissues during Surgery.

Electrical impedance spectroscopy (EIS) has been proposed as a promising noninvasive method to differentiate healthy thyroid from parathyroid tissues during thyroidectomy. However, previously reported similarities in the in vivo measured spectra of these tissues during a pilot study suggest that this separation may not be straightforward. We utilise computational modelling as a method to elucidate the distinguishing characteristics in the EIS signal and explore the features of the tissue that contribute to the observed electrical behaviour. Firstly, multiscale finite element models (or 'virtual tissue constructs') of thyroid and parathyroid tissues were developed and verified against in vivo tissue measurements. A global sensitivity analysis was performed to investigate the impact of physiological micro-, meso- and macroscale tissue morphological features of both tissue types on the computed macroscale EIS spectra and explore the separability of the two tissue types. Our results suggest that the presence of a surface fascia layer could obstruct tissue differentiation, but an analysis of the separability of simulated spectra without the surface fascia layer suggests that differentiation of the two tissue types should be possible if this layer is completely removed by the surgeon. Comprehensive in vivo measurements are required to fully determine the potential for EIS as a method in distinguishing between thyroid and parathyroid tissues.

Determining a musculoskeletal system's pre-stretched state using continuum-mechanical forward modelling and joint range optimization.

The subject-specific range of motion (RoM) of a musculoskeletal joint system is balanced by pre-tension levels of individual muscles, which affects their contraction capability. Such an inherent pre-tension or pre-stretch of muscles is not measureable with in vivo experiments. Using a 3D continuum mechanical forward simulation approach for motion analysis of the musculoskeletal system of the forearm with 3 flexor and 2 extensor muscles, we developed an optimization process to determine the muscle fibre pre-stretches for an initial arm position, which is given human dataset. We used RoM values of a healthy person to balance the motion in extension and flexion. The performed sensitivity study shows that the fibre pre-stretches of the m. brachialis, m. biceps brachii and m. triceps brachii with 91 % dominate the objective flexion ratio, while m. brachiradialis and m. anconeus amount 7.8 % and 1.2 % . Within the multi-dimensional space of the surrogate model, 3D sub-spaces of primary variables, namely the dominant muscles and the global objective, flexion ratio, exhibit a path of optimal solutions. Within this optimal path, the muscle fibre pre-stretch of two flexors demonstrate a negative correlation, while, in contrast, the primary extensor, m. triceps brachii correlates positively to each of the flexors. Comparing the global optimum with four other designs along the optimal path, we saw large deviations, e.g., up to 15 ∘ in motion and up to 40% in muscle force. This underlines the importance of accurate determination of fibre pre-stretch in muscles, especially, their role in pathological muscular disorders and surgical applications such as free muscle or tendon transfer.

Microcapsule Triggering Mechanics in Cementitious Materials: A Modelling and Machine Learning Approach.

Self-healing cementitious materials containing microcapsules filled with healing agents can autonomously seal cracks and restore structural integrity. However, optimising the microcapsule mechanical properties to survive concrete mixing whilst still rupturing at the cracked interface to release the healing agent remains challenging. This study develops an integrated numerical modelling and machine learning approach for tailoring acrylate-based microcapsules for triggering within cementitious matrices. Microfluidics is first utilised to produce microcapsules with systematically varied shell thickness, strength, and cement compatibility. The capsules are characterised and simulated using a continuum damage mechanics model that is able to simulate cracking. A parametric study investigates the key microcapsule and interfacial properties governing shell rupture versus matrix failure. The simulation results are used to train an artificial neural network to rapidly predict the triggering behaviour based on capsule properties. The machine learning model produces design curves relating the microcapsule strength, toughness, and interfacial bond to its propensity for fracture. By combining advanced simulations and data science, the framework connects tailored microcapsule properties to their intended performance in complex cementitious environments for more robust self-healing concrete systems.

Numerical model of pressure generated by elastic compression garments on a compressible human limb analogue.

OBJECTIVE: This study aimed to formulate a numerical approach (finite element modelling (FEM)) to calculate pressure values generated by compression garments on a compressible limb analogue, and to validate the numerical approach using experimental measurements. Existing models were also compared. METHOD: Experimentally measured pressure values and deformation caused by compression bands on a compressible human limb analogue were compared with values predicted using the Young-Laplace equation, a previously formulated analytical model and the FEM. RESULTS: The FEM provided greater accuracy in predicting the pressure generated by compression bands compared to existing models. The FEM also predicted deformation of the limb analogue with good agreement relative to experimental values. CONCLUSION: It was concluded that modelling the non-uniform manner in which the way a limb analogue is compressed should be incorporated into future modelling of the pressures generated by compression garments on a compressible limb analogue. DECLARATION OF INTEREST: The authors have no conflicts of interest to declare.

Numerical Modelling of the Thermoforming Behaviour of Thermoplastic Honeycomb Composite Sandwich Laminates.

Lightweight component design is effectively achievable through sandwich structures; many past research studies in the aerospace and racing sectors (since the 1920s) have proven it. To extend their application into the automotive and other transport industries, manufacturing cycle times must be reduced. This can be achieved by sandwich materials made of continuous fibre-reinforced thermoplastic (CFRTP) cover layers and thermoplastic honeycomb cores. To widen the application of flat thermoplastic-based sandwich panels into complex parts, a novel forming technology was developed by the Fraunhofer Institute of Microstructure of Materials and Systems (IMWS). Manufacturing defects like wrinkling and surface waviness should be minimised to achieve high reproducibility of the sandwich components. Studying different manufacturing parameters and their influence on the final part is complex and challenging to analyse through experiments, as it is time-consuming. Therefore, a finite element (FE) modelling approach is implemented to reduce such efforts. Initially, a thermoforming model is developed and validated with experimental results to check its reliability. Further, different simulations are performed to optimise the novel sandwich-forming process. In this study, a thermoplastic sandwich made of polypropylene (PP) honeycomb core and polypropylene glass fibre (PP/GF) cross-ply as cover layers was used, and its numerical model was executed in LS-DYNA software release R11.2.1.

Impact Behavior and Residual Strength of PEEK/CF-Laminated Composites with Various Stacking Sequences.

Carbon fiber-reinforced composites are popular due to their high strength and light weight; thus, the structures demonstrate high performance and specific strength. However, these composites are susceptible to impact damage. The objective of this research was to study the behavior of carbon fiber-reinforced laminates based on a polyetheretherketone (PEEK) matrix with six stacking sequences under static and impact loading. Four-point bending, short-beam bending, drop weight impact, and compression after impact tests were carried out. The results were complemented with digital shearography to estimate the damaged areas. Finite element modeling served to assess the failure mechanisms, such as fiber and matrix failure, in different layers due to tension of compression. Three behavior pattern of layups under drop-weight impact were found: (i)-energy redistribution due to mostly linear behavior (like a trampoline) and thus lower kinetic energy absorption for damage initiation, (ii)-moderate absorption of energy with initiation and propagation of concentrated damage with depressed redistribution of energy in the material, (iii)-moderate energy absorption with good redistribution due to initiation of small, dispersed damage. The results can be used to predict the mechanical behavior of composites with different stacking sequences in materials for proper structural design.

The impact of femoral head size on the wear evolution at contacting surfaces of total hip prostheses: A finite element analysis.

Total Hip Arthroplasty has been a revolutionary technique in restoring mobility to patients with damaged hip joints. The introduction of modular components of the hip prosthesis allowed for bespoke solutions based on the requirements of the patient. The femoral stem is designed with a conical trunnion to allow for assembly of different femoral head sizes based on surgical requirements. The femoral head diameters for a metal-on-polyethylene hip prosthesis have typically ranged between 22 mm and 36 mm and are typically manufactured using Cobalt-Chromium alloy. A smaller femoral head diameter is associated with lower wear of the polyethylene, however, there is a higher risk of dislocation. In this study, a finite element model of a standard commercial hip arthroplasty prosthesis was modelled with femoral head diameters ranging from 22 mm to 36 mm to investigate the wear evolution and material loss at both contacting surfaces (acetabular cup and femoral stem trunnion). The finite element model, coupled with a validated in-house wear algorithm modelled a human walking for 10 million steps. The results have shown that as the femoral head size increased, the amount of wear on all contacting surfaces increased. As the femoral head diameter increased from 22 mm to 36 mm, the highly cross-linked polyethylene (XLPE) volumetric wear increased by 61% from 98.6 mm3 to 159.5 mm3 while the femoral head taper surface volumetric wear increased by 21% from 4.18 mm3 to 4.95 mm3. This study has provided an insight into the amount of increased wear as the femoral head size increased which can highlight the life span of these prostheses in the human body.

The Effect of Lens Shape, Zonular Insertion and Finite Element Model on Simulated Shape Change of the Eye Lens.

The process of lens shape change in the eye to alter focussing (accommodation) is still not fully understood. Modelling approaches have been used to complement experimental findings in order to determine how constituents in the accommodative process influence the shape change of the lens. An unexplored factor in modelling is the role of the modelling software on the results of simulated shape change. Finite element models were constructed in both Abaqus and Ansys software using biological parameters from measurements of shape and refractive index of two 35-year-old lenses. The effect of zonular insertion on simulated shape change was tested on both 35-year-old lens models and with both types of software. Comparative analysis of shape change, optical power, and stress distributions showed that lens shape and zonular insertion positions affect the results of simulated shape change and that Abaqus and Ansys show differences in their respective models. The effect of the software package used needs to be taken into account when constructing finite element models and deriving conclusions.

Physical and 3D numerical modelling of reinforcements pullout test.

This paper reports results of laboratory and 3D numerical modeled pull-out tests with steel ladders and polymeric strip reinforcements. These types of reinforcement are commonly used in reinforced soil walls constructed with concrete facing elements. Laboratory pull-out tests are required to determine accurate and realistic pull-out strength values considering the interaction of specific reinforcement and backfill materials under different confining pressures (i.e., trying to simulate the different reinforcement layer arrangements and load conditions in actual reinforced soil walls). International design Codes for reinforced soil walls provide default values for pull-out strength. However, in many cases, default values are too conservative and/or are not strictly specified for particular reinforcement types. Pull-out tests can be difficult and expensive to perform, thus not being common nor worth for the vast majority of reinforced soil wall projects. Consequently, calibrated numerical models can be useful to predict pull-out response under site-specific conditions, and provide further understanding of the mechanisms involved in the soil-reinforcement interaction. Details of the numerical approach, including relevant aspects of the soil-reinforcement interfaces, are described. Examples of calibrated numerical predictions for pull-out loads, displacements, and soil-dilatancy effects are presented. The influence of reinforcement, soil and interface stiffnesses is shown. Numerical results provide useful insight for future modelling works of the complex interaction between type-specific backfill materials and reinforcement element, relevant for investigation and/or practical design of reinforced soil walls.

Experimental and numerical datasets for benchmark tests on high strength steel I-section frames.

This data article presents experimental and numerical datasets for eight fixed-base, single storey, unbraced high strength steel welded I-section frames subjected to in-plane horizontal and vertical loading. A detailed description of the full-scale frame testing programme is provided in the related research article 'Benchmark tests on high strength steel frames'. The experimental dataset can be used to steer future research in full-scale structural testing and provide benchmark results that are suitable for the validation of finite element models and the development of system-level design approaches. In addition to the benchmark experimental frame data, all necessary details and data for shell finite element (FE) model validation using geometrically and materially nonlinear analysis (GMNIA) is presented. The general purpose FE software Abaqus was used. The dataset can be used as an illustrative example of GMNIA validation in accordance with EN 1993-1-14 with all relevant data for reproducibility provided.

Local buckling of cold-formed steel trapezoidal sheets: Data for finite element model validation.

Data is provided from a validation example for a finite element model of a cold-formed steel trapezoidal sheet. The sheet is subjected to bending, failing due to local buckling. The numerical model and the validation procedure are carried out according to the new Eurocode 3 prEN1993-1-14 Design assisted by finite element analysis. Detailed information concerning all aspects needed to reproduce the example is included: (i) the nominal and measured values of the sheet geometry; (ii) the measured material properties of the steel; (iii) the test setup of the validation experiments; (iv) the experimental results; (v) a complete description of the finite element model and solution procedure; and (vi) the finite element results. Additionally, data related to sensitivity studies on the numerical model is also presented, including the effect of the model domain, meshing, and imperfections (shape, magnitude, direction and combinations). Overall, the article aims to provide data and guidance to designers and researchers validating similar numerical models.

3D-DXA Based Finite Element Modelling for Femur Strength Prediction: Evaluation Against QCT.

Osteoporosis is characterised by the loss of bone density resulting in an increased risk of fragility fractures. The clinical gold standard for diagnosing osteoporosis is based on the areal bone mineral density (aBMD) used as a surrogate for bone strength, in combination with clinical risk factors. Finite element (FE) analyses based on quantitative computed tomography (QCT) have been shown to estimate bone strength better than aBMD. However, their application in the osteoporosis clinics is limited due to exposure of patients to increased X-rays radiation dose. Statistical modelling methods (3D-DXA) enabling the estimation of 3D femur shape and volumetric bone density from dual energy X-ray absorptiometry (DXA) scan have been shown to improve osteoporosis management. The current study used 3D-DXA based FE analyses to estimate femur strength from the routine clinical DXA scans and compared its results against 151 QCT based FE analyses, in a clinical cohort of 157 subjects. The linear regression between the femur strength predicted by QCT-FE and 3D-DXA-FE models correlated highly (coefficient of determination R2 = 0.86) with a root mean square error (RMSE) of 397 N. In conclusion, the current study presented a 3D-DXA-FE modelling tool providing accurate femur strength estimates noninvasively, compared to QCT-FE models.

Data on the axial response of steel tubes infilled with rubberised alkali-activated concrete.

The presented data cover experimental and numerical axial load-shortening results of steel tubes infilled with rubberised alkali-activated concrete. The experimental data are obtained from 36 concrete filled steel tube specimens with circular and square cross-sections, length-to-diameter/width ratios of 2 and 4, and three different rubber contents in the concrete infill. The data from the numerical assessment cover the axial load-shortening response of over 300 finite element models. These cover a wide range of concrete infill strengths and rubber contents, steel tube grades, specimen widths, and steel tube wall thicknesses. Detailed descriptions of the material and methods, experimental testing, and numerical modelling procedures are also provided. The data reported herein supports the discussion in the research article "Axial compressive behaviour of composite steel elements incorporating rubberised alkali-activated concrete," and in the case of the numerical parametric assessment, give for the first time the full axial load-shortening response of all the models considered.

A macro-micro FE and ANN framework to assess site-specific bone ingrowth around the porous beaded-coated implant: an example with BOX® tibial implant for total ankle replacement.

The use of mechanoregulatory schemes based on finite element (FE) analysis for the evaluation of bone ingrowth around porous surfaces is a viable approach but requires significant computational time and effort. The aim of this study is to develop a combined macro-micro FE and artificial neural network (ANN) framework for rapid and accurate prediction of the site-specific bone ingrowth around the porous beaded-coated tibial implant for total ankle replacement (TAR). A macroscale FE model of the implanted tibia was developed based on CT data. Subsequently, a microscale FE model of the implant-bone interface was created for performing bone ingrowth simulations using mechanoregulatory algorithms. An ANN was trained for rapid and accurate prediction of bone ingrowth. The results predicted by ANN are well comparable to FE-predicted results. Predicted site-specific bone ingrowth using ANN around the implant ranges from 43.04 to 98.24%, with a mean bone ingrowth of around 74.24%. Results suggested that the central region exhibited the highest bone ingrowth, which is also well corroborated with the recent explanted study on BOX®. The proposed methodology has the potential to simulate bone ingrowth rapidly and effectively at any given site over any implant surface.

Dataset on sinusoidally stiffened 3D printed steel plated structures.

The paper reports a series of experimental and numerical data of destructive stub column tests on additively manufactured steel parts stiffened by surface sinusoidal wave patterns. The specimens were made in 316L stainless steel and manufactured by selective laser melting (SLM). The experimental tests covered five tensile coupon tests, fourteen square hollow section (SHS) stub column tests and measurements of geometric imperfections of the stub columns. Numerical models incorporating the measured material and geometric properties were developed and analysed via GMNIA approach. The validity of the numerical models is demonstrated by their accurate replications of the load-end shortening responses of the tested specimens. The reported dataset will contribute to the stability design and characterisation of thin-walled steel plated structures with advanced stiffening patterns.

Pre-intervention myocardial stress is a good predictor of aortic valvoluplasty outcome for fetal critical aortic stenosis and evolving HLHS.

Fetal critical aortic stenosis with evolving hypoplastic left heart syndrome (CAS-eHLHS) causes biomechanical and functional aberrations, leading to a high risk of progression to hypoplastic left heart syndrome (HLHS) at birth. Fetal aortic valvuloplasty (FAV) can resolve outflow obstruction and may reduce progression risk. However, it is currently difficult to accurately predict which patients will respond to the intervention and become functionally biventricular (BV) at birth, as opposed to becoming functionally univentricular (UV). This prediction is important for patient selection, parental counselling, and surgical planning. Therefore, we investigated whether biomechanics parameters from pre-FAV image-based computations could robustly distinguish between CAS-eHLHS cases with BV or UV outcomes in a retrospective cohort. To do so we performed image-based finite element biomechanics modelling of nine CAS-eHLHS cases undergoing intervention and six healthy fetal control hearts, and found that a biomechanical parameter, peak systolic myofibre stress, showed a uniquely large difference between BV and UV cases, which had a larger magnitude effect than echocardiography parameters. A simplified equation was derived for quick and easy estimation of myofibre stress from echo measurements via principal component analysis. When tested on a retrospective cohort of 37 CAS-eHLHS cases, the parameter outperformed other parameters in predicting UV versus BV outcomes, and thus has a high potential of improving outcome predictions, if incorporated into patient selection procedures. Physiologically, high myocardial stresses likely indicate a healthier myocardium that can withstand high stresses and resist pathological remodelling, which can explain why it is a good predictor of BV outcomes. KEY POINTS: Predicting the morphological birth outcomes (univentricular versus biventricular) of fetal aortic valvuloplasty for fetal aortic stenosis with evolving HLHS is important for accurate patient selection, parental counselling and management decisions. Computational simulations show that a biomechanics parameter, pre-intervention peak systolic myofibre stress, is uniquely robust in distinguishing between such outcomes, outperforming all echo parameters. An empirical equation was developed to quickly compute peak systolic myofibre stress from routine echo measurements and was the best predictor of outcomes among a wide range of parameters tested.

A novel computational workflow to holistically assess total knee arthroplasty biomechanics identifies subject-specific effects of joint mechanics on implant fixation.

Computational studies of total knee arthroplasty (TKA) often focus on either joint mechanics (kinematics and forces) or implant fixation mechanics. However, such disconnect between joint and fixation mechanics hinders our understanding of overall TKA biomechanical function by preventing identification of key relationships between these two levels of TKA mechanics. We developed a computational workflow to holistically assess TKA biomechanics by integrating musculoskeletal and finite element (FE) models. For our initial study using the workflow, we investigated how tibiofemoral contact mechanics affected the risk of failure due to debonding at the implant-cement interface using the four available subjects from the Grand Challenge Competitions to Predict In Vivo Knee Loads. We used a musculoskeletal model with a 12 degrees-of-freedom knee joint to simulate the stance phase of gait for each subject. The computed tibiofemoral joint forces at each node in contact were direct inputs to FE simulations of the same subjects. We found that the peak risk of failure did not coincide with the peak joint forces or the extreme tibiofemoral contact positions. Moreover, despite the consistency of joint forces across subjects, we observed important variability in the profile of the risk of failure during gait. Thus, by a combined evaluation of the joint and implant fixation mechanics of TKA, we could identify subject-specific effects of joint kinematics and forces on implant fixation that would otherwise have gone unnoticed. We intend to apply our workflow to evaluate the impact of implant alignment and design on TKA biomechanics.